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We study the conditions for the formation and propagation of Korteweg-de Vries (KdV) 

solitons in nuclear matter. In a previous work we have derived a KdV equation from Euler 

and continuity equations in non-relativistic hydrodynamics. In the present contribution 

we extend our formalism to relativistic fluids. We present results for a given equation of 

Cji . state, which is based on quantum hadrodynamics (QHD). 
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1. Introduction 



^ I Long ago |lj it was suggested that Korteweg - de Vries solitons might be formed in the 

•jH ■ nuclear medium. In a previous work [^ we have updated the early works on the subject 
1^ ■ introducing a realistic equation of state (EOS) for nuclear matter. We have found that 
lO ■ these solitary waves can indeed exist in the nuclear medium, provided that derivative 
■^ ! couplings between the nucleon and the vector field are included. These couplings lead to 
an energy density which depends on the Laplacian of the baryon density. For this class 
of equations of state, which is quite general as pointed out in [3], perturbations on the 
nuclear density can propagate as a pulse without dissipation. 

During the analysis of several realistic nuclear equations of state, we realized that, very 
often the speed of sound c^ is in the range 0.15 — 0.25. Compared to the speed of light 
these values are not large but not very small either. This suggests that, even for slowly 
moving nuclear matter, relativistic effects might be sizeable. This concern justifies the 
extension of the formalism presented in P] . 

2. Hydrodynamics 

Euler equation and the continuity equation form the basis of hydrodynamics. In the 
non-relativistic regime and for a perfect fluid they are [2]: 

|\(„-.V).-^-(±)v. (1) 

^ + V.(pbS)=0 (2) 

where pb, M, h and v are the baryon density, the nucleon mass, the enthalpy per nucleon 
and the fluid velocity respectively. In the relativistic case they are [^f5] : 

Ov -> \ f -> 0'd\ 
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^(Pi?7) + V-(PB7^?) = (4) 

where 7, e and p are the usual Lorentz factor (7 = (1 — v'^)~^^'^), energy density and 
pressure respectively. We have deliberately written the above equations in a non-covariant 
way to make the comparison between the non-relativistic and relativistic cases easier. 
Using the definition of enthalpy per nucleon p] for a perfect fiuid we find that dp = psdh. 
Therefore, using the Gibbs relation at zero temperature e + p = fis Pb '^^ can rewrite ([3]) 
as: 

where ^b is the baryochemical potential. Since the enthalpy per nucleon may also be 
written as [2H7l 



k=^ (6) 

OpB 

it becomes clear that the "force" on the right hand side of the Euler equations ([1]) and 
([5]) will be ultimately determined by the equation of state, i.e., by the function e{pB)- 

3. KdV equation and the nuclear equation of state 

Equations (IT]) and ([5]) contain the gradient of the derivative of the energy density. If e 
contains a Laplacian of pe, i-e., £ oc ... + ...V'^pb + ..., then ([T]) and (jS]) will have a cubic 
derivative with respect to the space coordinate which will give rise to the Korteweg-de 
Vries equation for the baryon density. The most popular relativistic mean field models 
do not have higher derivative terms and, even if they have at the start, these terms are 
usually neglected during the calculations. 

In [2] we have added a new derivative term to the usual non-linear QHD [S], given by 

Cm ^ -^.^t^id^d-'V.h'i^ (7) 

rriy 

where, as usual, the degrees of freedom are the baryon field "0, the neutral scalar meson 

field (f) and the neutral vector meson field V^, with the respective couplings and masses. 

The new term is designed to be small in comparison with the main baryon - vector meson 

interaction term g^Tp'^^V^il). Folowing the standard steps of the mean field formalism we 

arrive at the following expression for the energy density: 
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+ —(M*-Mf + ^pBV'pB (8) 



where r] is the baryon spin-isospin degeneracy factor, M* stands for the nucleon effective 
mass (given by M* = M — gs4>o) and the constants 6, c Qs and g^ are taken from [8]. 

Although Eq. ([8]) was obtained with the help of a specific Lagrangian taken from [8] 
and a prototype Laplacian interaction ([7]), the above form of the energy density follows 
quite naturally from an approach based on the density functional theory [5], regardless of 
the form of the underlying Lagrangian. Thus KdV solitons are a general consequence of 
many-body dynamics. 



4. KdV solitons 



We now repeat the steps developed in p!p] and introduce dimensionless variables for 
the baryon density and velocity: 



. PB 

P= — 

Po 



V = — 



We next define the "stretched coordinates" ^ and r as in pJffTO] : 
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where i? is a size scale and cr is a small (0 < a < 1) expansion parameter chosen to be 
[TO]: 
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where u is the propagation speed of the perturbation in question. We then expand ([9]) 
around the equilibrium values: 

p=l + api + a^p2 + ... (12) 

V = avi + cr^f2 + . . . (13) 

After the expansion above ([I]), ([2]), (jlj) and ([5]) will contain power series in a (in practice 
we go up to 0"^). Since the coefficients in these series are independent of each other we 
get a set of equations, which, when combined, lead to KdV equations for pi. In the 
non-relativistic case we have obtained [2]: 
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with the analytical solitonic solution: 
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gv V 2po 

where pi = api and u is the velocity of propagation of the perturbation. The solution 
above is a bump with width A given by: 
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Now, following the same sequence of steps, the combination of ([5]) and 
similar KdV equation for the relativistic case: 

dr ^ " ''^^ di \2MclmtR^) d^ 
with the solution given by: 

„ , , ?>{u — c 
pi[x,t) = : 



(16) 
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with the condition pB = M. 

As a consitency check we take the non-relativistic limit, which, in this case, means 
taking a small speed of sound c^ — )■ 0. In this limit (3 — c^^) = 3, flT7|) reduces to flT^ and 
f lTSj) coincides with f lTSj) . 



5. Conclusions 

The existence of KdV solitons in nuclear matter lias potential applications in nuclear 
physics at intermediate energies [T] and also possibly at high energies. The experimental 
measurements of jet quenching and related phenomena performed at RHIC [TT] offer an 
unique opportunity of studying supersonic motion in hot and dense hadronic matter. 
With this scenario in mind we gave the first step in the adaptation of the KdV soliton 
formalism to the new environment. We have extended the results of our previous work [2], 
showing that it is possible to obtain the KdV solitons in relativistic hydrodynamics with 
an appropriate EOS. Taking the non-relativistic limit {cj. — > 0) we were able to recover 
the previous results. 
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